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ABSTRACT 

Precise proper motion measurements (cr^ ~ 10/iasyr _1 ) of the recently discovered hyper- velocity 
star (HVS) SDSS J090745. 0+024507 would yield significant constraints on the axis ratios and ori- 
entation of a triaxial model for the Galactic halo. Triaxiality of dark matter halos is predicted by 
Cold Dark Matter models of galaxy formation and may be used to probe the nature of dark matter. 
However, unless the distance to this star is determined to better than 10%, these constraints suffer 
from one-dimensional degeneracies, which we quantify. We show how proper motion measurements 
of several HVSs could simultaneously resolve the distance degeneracies of all such stars and produce 
a detailed picture of the triaxial halo. Additional HVSs may be found from radial velocity surveys or 
from parallax/proper- motion data derived from GAIA. High-precision proper-motion measurements 
of these stars using the Space Interferometry Mission (SIM PlanetQuest) would substantially tighten 
the constraints they yield on the Galactic potential. 

Subject headings: cosmology: theory — dark matter : halos: structure — galaxies: formation 



1. INTRODUCTION 

In the course of a spectroscopic survey of candidate 
faint blue horizontal -branch stars found in the Sloan 
Digital Sky Survey, iBrown et alJ l)2005l) have discov- 
ered a star with the heliocentric radial velocity v h s = 
+853 ± 12 km s~ x at the Galactic coordinates b = 
31.3319°, I = 227.3353°. The velocity of this star, SDSS 
J090745. 0+024507 (hereafter called "the hyper-velocity 
star", or HVS), is more than twice the escape speed 
from the Galaxy at the star's present location. Cor- 
rected for the solar motion relative to the local standard 
of rest (LSR; iDehnen fe Binnevlll99^1 . the line-of-sight 
(los) LSR velocity is vi os — 848 km s _1 . The distance 
of the HVS from Earth is presently uncertain due to an 
ambiguity in its spectral classification. If it is a main 
sequence B9 star, its heliocentric distance is d « 70 kpc; 
if it is a blue horizontal branch star, then d ~ 40 kpc. 
A higher S/N spectrum from a large ground-based tele- 
scope should deter mine the spec t ral ty pe of this star. 

As discussed bv IBrown et all l)2005j) . the velocity of 
the HVS greatly exceeds that plausible for a runaway star 
ejected from a binary in which one component has under- 
gone a supernova explosion. The only known mechanism 
for a star to obtain such an extreme velocity is ejection 
from the deep potential of the massive black hole at the 
Galactic center, as a result of scattering with another star 
or tid al breakup of a binary l)Hillslll983 lYu fc Tremaind 
2003). Only extremely close to the massive black hole, 
at r < 0.01 pc, can stars attain the required speeds 
v w (2GM bh /r) 1 / 2 > 1000 km s" 1 . 

A precise measurement of the three-dimensional mo- 
tion of this HVS can probe the shape of the mass distri- 
bution of the Galactic halo in a way that is independent 
of any other technique attempted so far. The expected 
trajectory of the HVS in the Galaxy is almost a straight 
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line, owing to its extremely high velocity. However, as- 
suming its origin at the Galactic center, the direction of 
the HVS's present velocity should deviate slightly from 
being precisely radial due to departures from spherical 
symmetry of the Galactic potential. 

The main sources of such asymmetry are the flattened 
disk and the (possibly) triaxial dark matter halo. Sev- 
eral estimates of the halo shape based on observations of 
tidal debris associated with the Sagittarius d warf galaxy 
indicate that it is close to spherical (e.g. Ilbata et alJ 
l200lt iMaiewski etafl 1200! I.Iohnston et alJ 120051) . ah 
though Helmi (2004a) argues that minor-to-major axis 



ratios as low as 0.6 for isodensity contours cannot be 
ruled out for a prolate halo oriented with its major axis 
along the rotation axis of the disk. While cosmological 
N-body simulations based on Cold Dark Matter models 
typically produce prolate ha los with density axis ratios 
in th e range 0.5 - 0.8 (e.g. I.Iing fc Sutol l2002t iBullockl 
2002), gasdynamics simulations indicate that the effects 
of gas c ooling and dissipation t end to make the halos 
rounder ( Kaz antzidis et aT]l2004|) . Whether the predic- 
tions of gasdynamics simulations agree with the observa- 
tions of tidal streams still remains to be seen. 

Measuring the proper motions of the HVS and recon- 
structing its three-dimensional velocity will provide use- 
ful constraints on the shape and orientation of the Galac- 
tic dark matter halo, as we discuss below. Independent 
observational constraints of halo shapes are important 
for testing structure formation models as well as probing 
the nature of dark matter. 

2. ORBITS OF THE HYPER- VELOCITY STAR 

In order to evaluate the deviation of the HVS orbit 
from a straight line, we have calculated a family of orbits 
of the HVS consistent with its position in the sky (Galac- 
tic coordinates I and b), the assumed distance from Earth 
(d), and the observed line-of-sight velocity, v\ os . All or- 
bits start at the origin r = with some initial ejection 
velocity V{, which is then adjusted until the orbit repro- 
duces the four observables. 

We approximate the Galactic potential by the sum of 
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three components: spherical bulge, axisymmetric disk, 
and triaxial dark matter halo, with parameters that are 
consistent wi t h the current mass model of the Galaxy by 
Klvpin et al.l (12 002) . The bulge potential is given by the 



Hernauistl l|1990l ) model: 
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with mass Mf, = 10 M Q and core radius <z& = 0.6 
kpc. The disk p o tentia l is given by the analytical 
iMivamoto fc Nagail ()1975|) model: 



GM d 



(2) 



with mass Md = 4 x 10 M©, scale length a<j = 5 kpc, 
and scale height bd = 0.3 kpc. The halo potential is 
derived from a gen eralized triaxial den sity distribution 
of the NFW model ijNavarro et al.|[l99^) . which provides 
a good fit to halo profiles foun d in cosmological N-body 
simulations l|Jing fe Sutoll2002j) : 



Ph(f) 



(3) 



with mass Mh = 10 12 M and scale radius r s — 20 kpc. 
We generalize the halo profile as a function of the triaxial 
radius r to allow for three independent scale lengths: 



r 2 



{x/ qi f + {y/q 2 f + (z/q 3 f 
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We explore the full range of the axis ratios consistent 
with current observational constraints, 0.5 < qt < 1, 
to investigate the maximum effect of halo triaxiality on 
the expected proper motions of the HVS. Specifically, 
we calculate three sets of models with the halo major 
axis aligned with each of the three coordinate axes while 
varying the other two axis ratios from 0.5 to 1 with a 
step 0.05. For simplicity we assume that the axis ra- 
tios are constant as a function of radius, an assump- 
tion that is supported by the results of co smological N- 
body simulations (e.g. Ding fc Sutol l2002j) . In practice, 
the motion of the HVS probes the halo shape in the 
range of radii between 20 and 80 kpc (see Fig. [5J . We 
also do not include the effect of adiabatic contraction 
of the h alo in response to the central concentration of 
barvons dBlumenthal et~aDll986t iRvden fc Omul 119871 
Gncd in et al.ll2004l) since this effect is likely to be small 
at large radii, r > r s . 

The triaxial halo potential is obtained by numerical 
integration over thin tria xial shells, or homeoids (e.g. 
iBinnev fc Tremainelll987l Chapter 2.3). Interior to the 
shell the potential is constant, while isopotential surfaces 
exterior to the shell arc confocal ellipsoids labeled by the 
parameter r: 



const 



(5) 



The exterior potential due to each such shell m is pro- 
portional to ip(oo) — -0(m[r]), where 
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Fig. 1. — Deceleration of the HVS as a function of distance to the 
Galactic center. Three velocity components in the Galactocentric 
frame are shown for the cases of a spherical halo (dashed lines) and 
a triaxial halo (solid lines). The axis ratios of the triaxial model 
are qi/qz = 0.9, Q2/93 = 0.7, and the assumed current distance 
from Earth is d = 70 kpc. The bottom right panel shows the 
line-of-sight velocity as seen from Earth and the absolute value of 
the transverse velocity. The observed los velocity is shown by a 
horizontal dotted line. 



and the second equality is for the density profile given 
by equation ©. The total potential is the integral over 
all shells, 



®h(x,y,z) 



GM h g2g3 
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dr 
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We have tabulated the integral in equation Q, as well 
as the corresponding integrals for the force components, 
on a three-dimensional grid of coordinates (x, y, z) and 
interpolated the tables for orbit calculation. 

With our definition of the triaxial profile (eq. 0) the 
halo mass enclosed within a given spherical radius i? v i r 
depends on the axis ratios. The more triaxial halos have 
less enclosed mass. On the other hand, the measured 
mass of the Galaxy, inferred from the radial velocities 
of distant satellites and globular clusters, is in principle 
independent of halo triaxiality. In order to fix this dis- 
crepancy, we have renormalized the mass Mh for triaxial 
halo cases such that th e mass within a sp herical virial 
radius R V1Y — 12 r s (see iKlvpin et a,l.ll2002^) is the same 
for all models and equals 10 12 M . 

Figure ^ illustrates the deceleration of the HVS on its 
way from the Galactic center for both a spherical halo 
model and a triaxial halo model. A typical ejection ve- 
locity is 900 km s _1 in the Galactocentric frame, which 
is reduced to 700 km s _1 at the present position. The 
line-of-sight direction varies rapidly at small distances 
r < 20 kpc, when the HVS is relatively close to Earth, 
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Fig. 2. — Angle between the radius vector and the velocity 
vector of the HVS in the Galactocentric frame vs the distance to 
the Galactic center, for the orbits shown on Fig. Solid line 
shows the apparent deflection of the orbit from a straight line due 
to the triaxial halo, without the flattened disk. Dashed line shows 
the contribution of the disk alone, in the case of a spherical halo. 
At large distances from the center, most of the deflection is due 
to the asphericity of the halo. As a check of numerical accuracy, 
dotted line shows that in a spherical bulge+halo potential without 
the disk, the inferred angle is close to zero as expected. 



which leads to large apparent variations of the transverse 
velocity. At larger distances, a significant part of Vt T rep- 
resents the reflex motion of the Sun around the Galactic 
center. 

For this and subsequent figures, we have chosen a fidu- 
cial triaxial halo model with the axis ratios qi/q^ = 0.9, 
Q2/Q3 = 0.7, and the major axis aligned with the Z- 
coordinate. Our choice is motivated by the orientation 
of satellite galaxies and their counterparts in cosmologi- 
cal si mulations fe.g. JZentner et al.l2005HLibeskind et alJ 
2005) , which indicates that the dark matter halo may be 
prolate and oriented perpendicular to the plane of the 
disk. The prolate Galactic halo based on the kinemat- 
ics of t he Sag ittarius dwarf debris has been claimed by 
iHelmil l)2004b|) but disputed bv Uohnston et alJ (|2005^ . 
Even the eviden ce based on the sate l lites is inconclusive 
at present (c.f., INavarro et alJl200l IBailin et~all I2005|) 
and our fiducial model should be considered as an exam- 
ple only. 

Figure [21 shows the deviation of the orbit from a 
straight line, as measured by the angle between the ra- 
dius vector and the velocity vector in the Galactocentric 
frame. Owing to the very large ejection velocity from the 
center, this deviation is always small, < 1%. Therefore, 
the expected transverse velocities should be of the order 
1% of the total velocity, or several km s" 1 . 

Figure |21 also shows that the asymmetry of the poten- 
tial due to the flattened disk causes a smaller deflection 
than that due to the triaxial halo, for stars at large dis- 



FlG. 3. — Difference between the velocities perpendicular to the 
line of sight (solid line) for the orbits in the triaxial and spherical 
halos shown on Fig. For comparison, dashed line shows the dif- 
ference between the transverse velocities in the cases of a spherical 
halo with and without the flattened disk. These differences in Vtr 
lead to the differences in the expected proper motions of the HVS 
for different shapes of the Galactic halo. 



tances. In the case of a spherical halo, the deflection 
contributed by the disk peaks at r « 10 kpc but quickly 
declines at larger distances where the disk density van- 
ishes and the direction of the orbit aligns with the ve- 
locity vector. Because of this geometric effect, the angle 
#disk falls inversely proportional to the distance at r > 40 
kpc. On the other hand, the halo density profile is still 
close to isothermal at these distances, and the potential 
quadrupole induced by halo triaxiality, d$> oc ph(r)r 2 , is 
a weak function of r. This quadrupole causes a continu- 
ous deflection of the orbit, Sv ~ 5<$>/v. Hence, measuring 
the deviation of the trajectory of the HVS from a radial 
direction is sensitive to the halo triaxiality and relatively 
insensitive to any uncertainties in the mass model of the 
disk. 

A measurement of the proper motion of the HVS would 
give the transverse velocity in the reference frame asso- 
ciated with Earth (or after appropriate corrections, with 
the LSR). In addition to the perpendicular component 
of the Galactocentric velocity v, the measured velocity 
would include the transverse component of the Solar mo- 
tion around the Galaxy, v & , as well as a correction due 
to the Earth being at a distance Rq — 8 kpc from the 
Galactic center. At large distances, d 3> Rq, where the 
angle 8 <C 1 , the observed transverse velocity is 



^tr.obs 



vO 



Rq X ej_ 



HrQ, 



(8) 



where 8 is the vector of length 6 perpendicular to the 
line of sight, Rq is the vector of length i? Q towards the 
Galactic center, and e± is the unit vector perpendicular 
to both Rq and the line of sight. The last two compo- 
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Fig. 4. — Expected proper motions of the HVS under a range 
of different assumptions about the shape and orientation of the 
Galactic dark matter halo. The two components of the proper 
motion are in the directions of increasing Galactic coordinates, I 
and b. The family of models with the halo major axis along the 
Galactic X-coordinate is shown by triangles, along Y-coordinate 
by open squares, and along Z-coordinate by filled circles. The 
upper and lower groups of the proper motions are for two assumed 
distances of the HVS of 70 kpc and 40 kpc, respectively. The 
solid line shows the predicted proper motions for a HVS distance 
from Earth varying from 30 (bottom) to 140 kpc (top), assuming 
a spherical dark matter halo. 



nents depend only on the distance and angular coordi- 
nates of the HVS but not on the shape of the halo and 
the disk. For a given position of the HVS, the differ- 
ences in the expected proper motions are only due to the 
deflections by the halo and the disk. 

Figure shows the differences of the velocities perpen- 
dicular to the line-of-sight for the triaxial and spherical 
halos. For the chosen halo model, Au 4r « 6 km s _1 . 
This velocity difference is again dominated by the tri- 
axiality of the halo, as is the case with the deflection 
angles. At r > 1 kpc, the contribution of the flattened 
disk decreases as A«tr,disk <x d . The halo contribu- 
tion also decreases at large radii but only because of the 
smaller projection of the Solar reflex motion along the 
orbit (third term in eq. [H])- By construction of the or- 
bits that reproduce the observed positions of the HVS, 
at its current distance Avt r is independent of the solar 
motion. 

3. RECONSTRUCTION OF THE HALO SHAPE 

Placing constraints on the shape and orientation of 
the Galactic halo is important for testing current mod- 
els of galaxy formation. Current ground-based measure- 
ments of the proper motion of the HVS are inconclusive 
and consistent wit h zero within large uncertainties (see 
iBrown et alJl2005|) . However, the exquisite resolution of 
the Hubble Space Telescope (HST) would allow one to 



measure the proper motion of the HVS, or place useful 
limits on its value, within just a few years. 

The predicted proper motions in Galactic angular co- 
ordinates are plotted in Figure^J The velocity transverse 
to the line of sight is projected on two axes in the direc- 
tions of increasing Galactic coordinates I and b. The two 
groups of points represent sets of triaxial halo models at 
two discrete distances, while the solid curve represents a 
spherical model at a continuum of distances. The abso- 
lute values of the proper motions vary from 400/iasyr 
to 900/iasyr~ 1 , depending on which distance is chosen 
for the HVS (70 and 40 kpc, respectively). 

Note that the loci of the two distributions of the proper 
motions for two assumed distances are well separated 
from each other, and therefore the true distance of the 
HVS can be cleanly determined from a proper-motion 
measurement accurate to <7 M ~ 100/iasyr -1 . 

With an accuracy of <7 M ~ 20/iasyr -1 we can place 
interesting constraints on the orientation of the Galactic 
halo. In particular, if the major axis of the halo is aligned 
with the direction to the Galactic center (I = 0), the 
predicted components m all lie below m < 280^asyr~ 1 . 
If the major axis is in the direction of solar rotation, fj,i > 
SOO^asyr" 1 . If the halo is prolate with the major axis 
aligned with the disk rotation axis, then the predicted 
fib component of the proper motion is well constrained 
to be either -320 ± 20 or -700 ± 40^asyr~ 1 , for the 
assumed distances of 70 and 40 kpc, respectively. These 
estimates are valid under our assumption that the halo 
axes are aligned with the Galactic coordinate axes, i.e., 
X along I = 0, Y along / = 90°, and Z perpendicular to 
the disk plane. 

With an accuracy of <7 M = 10/iasyr" 1 we can attempt 
a reconstruction of the halo axis ratios. Figure shows 
an example of such a future reconstruction given a hypo- 
thetical measurement of the proper motions of the HVS, 
fi° hs and /4J bs . The three sets of models (each with the 
major axis aligned with one of the coordinate axes) are 
combined in a single plot of the axis ratios (72 / q\ = Y/X 
versus (73/(71 = Z/X. The lower left part of the plot cor- 
responds to the case of the major axis along the Galactic 
X-coordinate (qi = 1), whereas the lower right and up- 
per left parts correspond to the major axis along the 
Y-coordinate ((72 = 1) and the Z-coordinate ((73 = 1), 
respectively. In this plot we calculate 
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assuming that the true distance of the HVS is unknown 
(a e i — > 00) and search for all predicted proper motions 
[iii [ib a t au possible distances between 20 kpc and 100 
kpc that minimize \ 2 . The assumed observed proper 
motions, however, correspond to the fiducial model with 
?i/?3 = 0.9, (72/(73 = 0.7, and the distance of 70 kpc. 

According to FigureEl a certain range of the axis ratios 
can be excluded based on the proper motion accuracy of 
= lO/iasyr -1 . However, there are strong triaxiality- 
distance degeneracies. The best-fitting axis ratios form 
a band stretching across the diagram because different 
halo shapes with a somewhat different target distance 
can produce similar proper motions. This effect is evi- 
dent in Figure The inferred distances for our fiducial 
model range from 67 to 78 kpc. The allowed contours 
would be reduced if the actual distance of the HVS were 
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Fig. 5. — Reconstruction of the halo shape from the (assumed) 
measured proper motions of the HVS. All calculated triaxial models 
are marked by filled dots corresponding to their axis ratios. Lines 
represent la, 2a, and 3a contours, assuming measurement errors 
of 10/iasyr — 1 . The assumed proper motions correspond to the 
triaxial model with the axis ratios qi/qs = 0.9, 52/93 = 0.7, which 
are marked by a star, and the distance of the HVS of d = 70 kpc. 



known, but it needs to be measured to better than 10% 
to make a significant difference. Also, if the true distance 
of the HVS is smaller than the 70 kpc assumed in our 
fiducial model, the constraints become tighter. 

4. TWO, THREE, ... MANY HYPER VELOCITY STARS 
4.1. Two Hyper Velocity Stars 

Measurements of the proper motions of the single 
known HVS would already provide interesting con- 
straints on the shape and orientation of the Galactic halo. 
Yet, the allowed range of axis ratios is fairly broad, even 
when the measurement errors are as small as 10 /lasyr -1 . 
The constraints could become tighter if another HVS 
were discovered. 

For example, suppose that another hypervelocity star, 
HVS2, were discovered at the same angle 6 2 = 31° above 
the Galactic plane but offset by 45° in longitude, I2 = 
I — 45° = 182°, and at a hypothetical distance of 50 kpc. 
Figure El shows the superposition of the contours of the 
original HVS and the new HVS2. The x 2 contours for 
the two stars intersect almost perpendicularly to each 
other and significantly reduce the acceptable range of 
halo shapes. At the la level the axis ratios would be 
constrained to better than 20%. 

Do all directions of the HVS2 on the sky provide 
equally useful constraints on the halo axis ratios? Figure 
□ shows the contours for the HVS2 offset by 180° in lon- 
gitude (l 2 = 47°) and 30° in latitude (b 2 = 61°). These 
contours also intersect the contours for the original HVS 
and the resulting constraints are also tight. However, 
in the case (not shown) that the two stars were at the 



(b,l,d) = (31° 227° 70 kpc) and (31° 182° 50 kpc) 




axis ratio Y/X 

Fig. 6. — Combined reconstructed contours for two HVS, as- 
suming another star were detected at a 45° offset in the Galactic I 
coordinate, at a distance of 50 kpc from Earth. 



(b,l,d) = (31° 227° 70 kpc) and (61° 47° 50 kpc) 




0.5 1 1.5 2 



axis ratio Y/X 

Fig. 7. — Same as Fig. [5] but assuming an offset of 180° in £ and 
30° in b. 



same longitude (I2 = 227°) but only offset in latitude 
(b-2 = 61°), the contours of the HVS2 are almost paral- 
lel to those of the HVS and provide no new information. 
Thus, searching prospective HVSs at different longitudes 
than the original HVS would yield the most interesting 
new constraints. 
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4.2. Three or more Hyper Velocity Stars 

Although there are physical motivations (albeit not 
conclusive) for the halo to be prolate with its major 
axis perpend i cular to the plane o f the Galactic disk (e.g. 
lHehnill2004allZentner et al.ll2005|) . the orientation of the 
other two axes in the plane of the disk is somewhat ar- 
bitrary. There is nothing special about the current loca- 
tion of the Sun, which defines the Galactic A-coordinatc. 
Therefore, the halo axes in the plane can be rotated with 
respect to the Galactic coordinates by any angle between 
and 90°. This means that even under the most opti- 
mistic assumptions, a minimum of three HVSs are re- 
quired to fix the Galactic potential. 

To illustrate the role of three HVSs, we have calculated 
another set of models, keeping fixed the ratio qi/qz = 0.9 
and varying two parameters: the ratio 92/93 between 0.5 
and 0.9, and the angle <fi between the halo axis q\ and 
the Galactic A-coordinate between 5° and 90°. Figure 
[S] shows the resulting contours for a fiducial model in 
which the qi axis is rotated by 60°. The rotation angle is 
constrained very well, within 5° — 10°, even if we assume 
larger measurement errors, = 30/xasyr -1 . 

We have also explored the configuration (not shown) 
in which the halo is oblate rather than prolate, with the 
minor axis along the rotation axis of the disk. In this 
configuration the angular momentum vectors of the halo 
and the disk are parallel to each other, as has often been 
assumed in the literature. We find that the contours 
in this model are qualitatively similar to those in Fig. 
[8] and in particular, angle <p can be constrained with a 
similar accuracy. 

If five HVSs were detectable, it would be possible to 
measure all five parameters of a more general triaxial 
halo. These would be the two axis ratios, and the three 
direction angles. Additional HVSs (beyond five) could 
be used either to improve the precision of the estimates 
of the first five parameters or to explore halo models of 
additional complexity. 

5. DISCUSSION 

We have shown that proper motion measurements 
of the recently discovered HVS, accurate to cr^ ~ 
lO^asyr -1 , would place significant constraints on the 
axis ratios and orientation of the Galactic halo. With 
only one HVS, such constraints are limited by the 
triaxiality-distance degeneracies. They can be improved 
only when the distance is determined to within a few 
kpc. However, the discoveries and proper motion mea- 
surements of several additional HVSs could simultane- 
ously resolve the distance degeneracies of all such stars 
and produce tighter limits on the triaxiality of the Galac- 
tic halo. We now discuss implications of the discovery of 
multiple HVSs and potential methods of finding more 
HVSs in future surveys. 

5.1. History of Ejection from the Galactic Center 

If several hyper velocity stars are found, it may be 
possible to reconstruct the history of their ejections 
from the Galactic center and therefore, constrain the 
rate of binary disrupt ion near the massive black hole. 
lYu fc Tremainei l)2003l) estimate the rate of ejection from 
the breakup of binary stars of ~ 10~ 5 yr _1 . If there is 
another massive black hole in the close vicinity of the 
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halo angle </> 

Fig. 8. — Same as Fig. but as a function of angle between 
the halo axis qi and the Galactic X-coordinate. The axis ratio 
9l/93 = 0-9 i s kept fixed in this plot. Assumed measurement errors 
are 30 /ias yr _1 . 



central black hole at the Galactic center, then stars can 
be ejected by three-body interactions with the binary 
system formed by the two black holes. lYu fc Tremainei 
( 2003) estimate that the rate of such events < 10 -4 yr -1 . 
The density of the HVSs in the halo can be calculated if 
HVSs are ejected from the Galactic center isotropically 
at a steady rate N = 10~ 5 AV5 yr -1 . The expected num- 
ber density in this case is n(r) « A_5 r~ 2 kpc -3 , where 
r is in kpc. 

There are several circumstances that could lead to 
clumping of HVS ejections in time, e.g., if bursts of 
star formation (lasting ~ 1 Myr) or infall of clus- 
ters containing intermediate -mass black holes (e.g., 
IHansen fc^ Milosavlicvic 200j|) led to a rapid increase in 
the ejection rate. If the global fit to multi-HVS data 
led to distance estimates accurate to 5 kpc, it would en- 
able one to date individual ejections with an accuracy of 
about 5 Myr, and therefore, to test whether the ejections 
were clumped in time. 

5.2. How to find HVS? 

There are basically two potential methods to find more 
HVSs: proper motions and radial velocities. For prac- 
tical reasons, these methods probe two different vol- 
umes, the (~ 10 kpc) "solar neighborhood" and the outer 
Galaxy, respectively. 

If a star ejected from the Galactic center is pass- 
ing within ~ 10 kpc of the Sun with a speed signifi- 
cantly greater than the escape speed, then its proper 
motion will be of order fj, ~ 1000kms _1 /10kpc ~ 
20masyr~ 1 . Given a proper motion survey with preci- 
sion ~ 4masyr~ 1 , this would be detectable at the 5ct 
level. Such precision is already available for the sev- 
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eral thousand sq uare degrees covered by the Sloan Dig- 
italSky Survey HConld KollTrieierll200i IMnnn et all 
2004). However, robust identification of high- velocity 
candidates requires not just reliable proper motions, 
but reliable distances as well. While horizontal branch 
(HB) stars are approximately standard candles whose 
distances can be reliably estimated from their colors and 
magnitudes, halo blue stragglers (BS) can masquerade 
as HB stars. A BS that was 2 mag dimmer than a HB 
star would have the same proper motion as a HB HVS 
but would actually be travelling 2.5 times slower. That 
is, an ordinary halo BS could easily be misinterpreted 
as a HB HVS. Similarly, because subdwarfs are often 2 
mag dimmer than main-sequence stars of the same color, 
ordinary subdwarfs could be mistaken for HVS main- 
sequence stars. To some extent, HVS candidates can be 
singled out because their trajectories "point back" to the 
Galactic center. However, since halo stars tend to be on 
radial orbits, this characteristic is not so distinguishing 
unless the proper motions are fairly precise. 

The situation will change radically when GAIA data 
become available, not only because of its An sky cover- 
age but also because it will provide parallaxes as well 
as proper motions. At present, these are expected to be 
ov ~ 20 ^as at V = 15.5, implying roughly 20% paral- 
lax errors for HB stars at 10 kpc. Dimmer stars would 
be similarly measurable but only within a much smaller 
volume. 

The proper-motion technique does not work at all for 
HVS stars like SDSS J090745.0+024507 that have al- 
ready left the Galaxy. The transverse velocities of these 
stars in the Galactic frame are only w tr ~ v(Ro/d), 
where v is their actual velocity and d is their distance 
(see eq. [HI). This is similar to the value for halo 
stars, making them not easily distinguishable. Rather, 
these distant HVS can be reliably identified using radial- 
velocity surveys of distant stars (as was done for SDSS 
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J090745.0+024507). The only drawback is that many 
stars with large photometric distances must be measured 
to find a small number of HVSs. Since the number den- 
sity of HVS falls as r~ 2 , while the density of contam- 
inating Galactic stars falls as r~ 3 5 , the contamination 
rate falls for more distant samples. However, more dis- 
tant stars require longer integration times. Moreover, 
once found, they require more accurate proper motions 
(on fainter stars!) to derive the same precision on the 
transverse velocity, and so to extract similar information 
about the Galactic potential. Thus, there do not appear 
to be any easy paths to finding HVSs. 

Finally, we remark that the simulations we have con- 
ducted have assumed proper-motion errors of ~ 10 — 
50/iasyr -1 , the upper end of which is roughly what is 
achievable using the HST Advanced Camera for Surveys 
with a 4-year baseline. It is open to question whether 
HST will even be operational in 4 years, let alone in the 
additional time required to first find the HVSs. Over 
the longer term, the Space Inter ferometry Mission (SIM 
PlanetQuest) expects to achieve c M ~ 4//asyr~ 1 for tar- 
geted stars with V < 20. For stars toward this faint limit, 
this very high precision can be achieved only by quite 
long integrations, typically several tens of hours. How- 
ever, the integration time falls inversely as the square 
of the desired precision, so that exposure times could 
be fine-tuned to the precision required for each specific 
HVS. 
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